The phase diagram of high-Tc's: Influence of anisotropy and disorder 
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We propose a phase diagram for the vortex structure of high temperature superconductors which 
incorporates the effects of anisotropy and disorder, ft is based on numerical simulations using the 
three-dimensional Josephson junction array model. We support the results with an estimation of 
the internal energy and configurational entropy of the system. Our results give a unified picture of 
the behavior of the vortex lattice, covering from the very anysotropic Bi2Sr2CaCu208 to the less 
anisotropic YBa2Cu307, and from the first order melting ocurring in clean samples to the continuous 
transitions observed in samples with defects. 
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I. INTRODUCTION 



The phase diagram of high temperature superconduc- 
tors in the mixed state has provided an astonishingly 
broad field to workers in -among other fields- many 
body problems, polymer physics, low dimensional sys- 
tems, critical phenomena, and statistical physics in gen- 
eral. The main reason of this situation is the great num- 
ber of parameters that define the behavior of the vor- 
tex structure. On the other hand, the same abundance 
of parameters defining the system turns difficult to find 
a unified description of all features observed in experi- 
ments. Some of the main parameters that define the be- 
havior of the vortex structure are the external magnetic 
field H, temperature T, anisotropy rj, and the disorder 
(which produces a non-homogeneous pinning potential 
for the vortices) that at this moment we loosely charac- 
terize by a parameter D. There are convincing explana- 
tions of the main characteristics of different sectors of this 
multi-dimensional phase diagram, such as the first order 
melting of the vortex-lattice in clean samples, the con- 
tinuous melting of a glassy phase in disordered samples, 
or the existence of two different superconducting tran- 
sitions (perpendicular and parallel to H) in some cases 
(for a review see Ref. 1). However a unified, consistent 
with experiments description of the problem, even in a 
qualitative level, is still lacking. 

In this paper we propose a qualitative H-T-D-rj phase 
diagram of high-T^s, that reproduces most of the avail- 
able experimental results. Our approach is twofold: 
we use numerical simulation on the three-dimensional 
Josephson junction array model to study the behavior 
of the system as a function of D and 77, and show that 
the dependence on H can be deduced from a rescaling of 
D and 77. The obtained phase diagram is then rederived 
using a phenomenological estimation of the free energy F 
of the system for different values of D and rj. This estima- 
tion relies on the existence of two characteristic lengths 
and ^ab parallel and perpendicular to the applied field H 
which are supposed to govern the behavior of the system. 



The minimizing of F with respect to and ^ab allows 
one to obtain the S,c{T) and S,ab{T) functions, which in 
turn are used to detect the superconducting transitions. 

The reminder of the paper is organized as follows. In 
the next section we present the results of the numerical 
simulations, and discuss the D-rj phase diagram emerg- 
ing from them. In Sec. Ill this phase diagram is quali- 
tatively re-obtained using a proposal for the free energy 
of the system. In Sec. IV we indicate that a change in 
the external magnetic field can be interpreted as a move- 
ments in the D-t] plane, and so the H-T phase diagram 
for samples with different D and 77 can be obtained from 
the results of the previous sections. We also compare our 
results with those found in experimental studies. Finally 
in Sec. V, we summarize and conclude. 



II. NUMERICAL SIMULATIONS 

A. The model 

Our numerical results are based on simulations per- 
formed on the three dimensional (3D) Josephson junc- 
tion array (JJA) model on a stacked triangular network. 
Each junction is modelled by an ideal Josephson junc- 
tion with critical current Ic shunted by a normal resis- 
tance R and its attached Johnson noise generator, which 
accounts for the effects of temperature. An external c 
directed magnetic field is included. The variables that 
characterize the model are the phase of a superconduct- 
ing order parameter defined on the nodes of the lattice. 
Vortices form in the system as singularities in the dis- 
tribution of these phases- The details of the model have 
been discussed elsewhere.B'cl The 3D JJA model has been 
previously used to show the ficat order melting of the 
vortex lattice in clean systems.aEl Both thermodynami- 
cal and transport signatures of this first order transition 
were obtained, in close relation to experimental results.u 
In addition, using the same model we have siown that 
disorder can destroy the first order transition^ 
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We study here the model further by systematically ex- 
ploring the case of anisotropic and disordered samples. 
We introduce anisotropy by reducing the critical current 
of the c axis directed junctions by a factor 77^, and at the 
same time increasing the c axis normal resistance by the 
same factor. Disorder is introduced by randomly varying 
the critical current of the junctions through the lattice. 
As vortices gain energy when close to a low critical cur- 
rent region, the effect of randomizing the critical currents 
is to provide a nearly random pinning potential for vor- 
tices. We characterize the disorder by a parameter D 
which is defined as D = {l'^'^'' - ''^) / (/^"^'^ + /™'") , 
where I'^^'^^ and are the maximum and minimum 
value of the critical current of the junctions through the 
sample. The probability distribution between J™'*'^ and 
If'"" is taken flat. 

We carried out simulations for H = 1/6 flux quanta 
per plaquette. This is the value used in Refs. 4 and 5. 
It produces a ground state (for a clean sample) which is 
commensurate with the subjacent triangular lattice, so 
no frustration effects are expected. Although the value 
1/6 is rather large and effects of the substrate may be 
observable, we expect the physics of the problem to be 
qualitatively well described. In particular, we assume 
that for a clean sample the first order melting observed 
in simulations is in jfapt-the counterpart of the experi- 
mental observations .Boiiil It would be interesting to per- 
form simulations at lower (commensurate) fields, such as 
1/14 or 1/36. However, the sample size needed to mini- 
mize size effects make the computing time be exceedingly 
large. 



B. Results 

All simulations presented here were performed for H = 
1/6, with boundary conditions as in Ref. 7, and for a 
sample of Lab x Lab x ic = 18 x 18 x 18 junctions. We 
characterize the superconducting transitions by measur- 
ing the resistivity of the sample when a small current 
(typically 1/100 of the mean critical current in the corre- 
sponding direction) is applied along the ab or c direction. 
We will observe two well different behaviors of the resis- 
tivities as a function of temperature: in some cases resis- 
tivities have a jump from zero to a finite value at a given 
temperature. This jump in the resistivity corresponds 
(see the discussions in Sec. IV) to a first order phase 
transition of the vortex lattice. In some other cases we 
will obtain that resistivities as a function of temperature 
smoothly depart from zero. We will refer to this behav- 
ior as a continuous transition. We do not claim at this 
point about whether these continuous transitions are or 
are not real phase transitions, they can be crossovers as 
well. In Sec. V we present a discussion on this point. 

We start by showing in Fig. |](a) and (b) (full symbols) 
the results for the ab plane and c axis resistivity of an 
isotropic sample as a function of temperature for differ- 



ent values of the disorder D. For D = we obtain a jump 
in the resistivities at a well defined temperature (the low 
temperature tails in pab are due to surface effects). This 
temperature is the (first order) melting temperature 
of the system. At T„i the superconducting coherence 
is lost discontinuously in all directions. In other words, 
vortices passes from a solid phase to an entangled liquid 
phase. This situation persists for low values of disor- 
der. When disorder increases further (D > 0.3 in our 
simulations), the first order transition is lost, and two 
continuous transitions at different temperatures Ti -for 
Pab- and Tp -for pc- are obtained. Superconducting co- 
herence along the ab plane is lost at Ti when increasing 
the temperature, but the system is still superconducting 
along the c direction. At the higher temperature Tp the 
parallel to field superconducting coherence is lost. For 
Ti < T < Tp the vortex structure corresponds to that of 
a disentangled vortex liquid. 
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FIG. 1. Resistivities pab, pc, helicity modulus Xab and Xc 
and percolation probability P as a function of temperature 
for ri — 1 and different values of the disorder D as indicated. 
Temperature is measured in units of the mean Josephson en- 
ergy of the in- plane junctions. Magnetic field is = 1/6 
quantum fluxes per plaquette. 



Thermodynamical measurements agree with this pic- 
ture. We calculated the helicity modulus parallel (xc) 
and perpendicular (xab) to the field. Helicity modulus 
measures the influence on the energy of the system of a 
twist in the boundary conditions, and has to be different 
from zero to indicate superconducting coherence.E3 Open 
symbols in Fig. ^a) and (b) show the values x that cor- 
respond to the resistivity curves. In the case of the first 
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order transition both Xc and Xab have an abrupt drop at 
the melting temperature. For the highly disordered case 
Xab has a (smoother) decrease at Ti, whereas Xc becomes 
nearly zero at Tp. The transition at Tp is a percolation 
phase transition of the vortex structure as discussed in 
Refs. 14 and 12, that can be characterized by a perco- 
lation probability P. The value of P is the probability 
that a vortex path traversing the sample perpendicularly 
to the applied field exists. This value is zero below the 
transition and one above, with a transition zone that be- 
comes narrower when Lab increases. The plot of P vs 
temperature is shown in Fig. |^(c). We see that pc is 
different from (equal to) zero if P is different from (equal 
to) zero. 
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FIG. 2. Same as Fig. 1 but for a sample with anisotropy 
rf = 2Q. The values of Xc for D = 0.3 and D = 0.5 are nearly 
zero within the considered temperature range. 



paths across the samples exist, however these path are 
not mobile (presumably because they are pinned to, the 
planes, as T < Ti) and dissipation is not observed.lla 

From the above discussion on Figs. |^ and ||, the fol- 
lowing picture emerges: at low disorder the vortex lattice 
melts through a first order phase transition. When dis- 
order increases this transition breaks into two continuous 
ones, where in-plane and inter-plane coherence are lost at 
different temperatures Ti and Tp. Whether Ti is larger or 
lower than Tp depends on the anisotropy of the system. 
For nearly isotropic samples Ti is lower than Tp, and a 
disentangled vortex liquid phase exists for Ti < T < Tp. 
For highly anisotropic samples Tp is lower than Ti al- 
though dissipation along the c axis is observed only for 
T>T,. 

In Fig. ^ we show the numerical D-i] phase diagram as 
obtained from simulations. Stars indicate points where 
the melting transition is first order. Full circles are 
points where Ti < Tp, and hollow circles are points where 
Tp < Ti. This diagram shows the three different zones 
discussed above. The continuous line is a sketch of the 
frontiers between the different zones. In the next section 
we show that the main characteristics of this phase dia- 
gram can be derived from a simplified description of the 
vortex structure. 
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In Fig. we give the corresponding results for 
anisotropy ry — 20. The first order melting transition 
occurring for low disorder is similar to that observed 
for T] = 1, i.e., the inter-plane coherence is lost at the 
same temperature than the in-plane coherence. How- 
ever, for higher values of disorder some differences occur: 
the inter-plane resistive transition occurs (to our numer- 
ical precision) at the same temperature as the in-plane 
transition Ti , but the percolation transition temperature 
Tp (which for low anisotropy coincides with the resistive 
inter-plane transition) occurs for lower temperatures, i.e., 
Tp < Ti. This indicates that for high anisotropy, there is 
a temperature range Tp < T < Ti for which percolation 



FIG. 3. Disorder-anisotropy phase diagram as obtained 
from the simulations. Stars indicate points where the melting 
transition is first order. Full circles are points where T <Tp, 
and hollow circles are points where Tp < Ti. The continuous 
line is a sketch of the frontiers between the different zones. 



III. A SIMPLIFIED MODEL 

We saw in the previous section that a single first order 
transition is observed at low disorder, whereas two con- 
tinuous ones are obtained in disordered samples. These 
features suggest that the system can be described as hav- 
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ing in general two different transitions, but that in certain 
cases tliey can merge onto a single first order transition 
due to some kind of 'interaction'. Here we show that this 
idea can be formulated more precisely by estimating the 
free energy of the system. 

A. Free energy functional 

To make this estimation we will consider first a sin- 
gle clean plane. We suppose that the thermodynamics 
of that plane can be phenomenologically described by a 
quantity ^ab, which is a correlation length: for distances 
shorter than the system has superconducting coher- 
ence, whereas this coherence is lost for distances larger 
than ^ab- The free energy functional J^"^^ of the plane 
has the form 

^2^ = £;2^(eab)-r52^(U) (1) 

(the thermodynamical free energy F is obtained by min- 
imizing with respect to ^ab)- For a system of com- 



pletely decoupled planes, we would have 

= L^E'^'iU) - TL^S^^iU) (2) 

On the other hand, if the coupling between planes is in- 
finite the vortices are rigid lines and we get 

^,'^0 = - {o^hf) Lc + ici?'^(U) - TS^^iU)- (3) 

Note that in this case the entropy term does not have 
the factor because giving the position of the vortices 
on one plane automatically determines the position of 
vortices in all other planes. The term — (a/ry^) Lc ac- 
counts for the energy gain due to the coupling of the 
planes, a being a numerical constant and ry the anisotropy 
parameter defined before. In an intermediate situation 
(0 < 7/ < oo) the system can be thought as formed by 
ic/^c layers (^c is a number that satisfies 1 < i^c < Lc). 
Within each layer the vortices are almost straight lines, 
whereas correlation is small between different layers. 
Within this picture the free energy of the system is 



- - (a/ry^) (L, - Lj^c) + ici?'^(U) - T (L./Cc) S^°{U). (4) 

In the first term, Lc — LcjS^c is the number of sites along the c-direction at which the system gains an energy 
a/rf. The length is a correlation length along the c-direction. For distances shorter than the system possesses 
superconducting coherence, whereas this coherence is lost at distances greater than ^c- 

The previous estimation of the free energy of the system is too crude. In particular, in the form given by Eq. (Q), 
it leads to some unphysical results. We must modify Eq. slightly in order to obtain the correct behavior in some 
limiting cases. However, we will keep a fundamental property of Eq. and make the guess that the entropy can be 
written for the real system as 

S^fi^S^'^iU), (5) 

i.e., as a product of independent functions of ^c and ^^fc, with /(^c) and S'^^{S_ab) two yet unknown functions. The 
basic assumption contained in (^) is the following: if the value of is kept fixed, then the system behaves as a two 
dimensional system with a renormalized temperature. Although this assumption cannot be fully justified a priori, it 
is a natural starting point, and gives sensible results as we will show soon. 

We still have to add a term to the free energy which accounts for the effect of impurities. Impurities decrease the 
energy of the system when vortices pin to them. If pinning is uncorrclated the energy gain due to pinning becomes 
lower when the vortex positions are more correlated. So we add a term to the free energy that increases the energy 
of the system when ^af, and increase. The most simple term of this type is of the form D^ab^c- Finally, the free 
energy functional of the system J-{^abT^c) (dropping an irrelevant constant term) is 

^ iU, = (a/v^) Lc/ic + LcE^'^iU) - Tf{^c)S^''{U) + DU^c (6) 

The true expressions for the functions /, E^^ and S^^ are difficult to establish. However, it is not our aim to give 
a complete quantitative description of the free energy of the system but only a qualitative description of the phase 
diagram. We will only ask the functions /, E^^ and S^^ to reproduce some known limiting cases: if (fab) is kept 
fix we expect the value of £,ab{T) {^c{T)) obtained by minimizing J-, to be smoothly dependent on tempecature. This 
corresponds to the absence of first order transitions in the dynamics of a single plane (a single vortex linc).llj However, 
when J- is minimized with respect to both ^c and fafc a discontinuity in S,c{T) and ^abiT) can appear, as we show 
below. 

Just to give an example we use for the function / the form / = 7 In (ic/Cc) + 1- The term 7 In (ic/Cc) (7 is a 
numerical constant) is proportional to the entropy of a single (isolated) vortex line. The constant added assures that 
the two-dimensional limiting case is reobtained when Lc = fc- In addition, we will take for the functions E^^ and 
S^^ the form E'^^ ~ Lab/^ab and S'^^ = 7ln {Lab/^ab) -I- 1 , in such a way that we obtain a form for the free energy 
fimctional that is (unphysically!) symmetric (for a/r/'^ — 1) between ab- and c-directions. We have tested other forms 



4 



of the functions /, E^^ and S^^ (giving the same limiting behavior discussed above) and found results qualitatively 
similar to those shown here. 

We arrive to our final working expression of the free energy functional (we will measure S^ab and in units of Lab 
and Lc, respectively, take a = 1 by rescaling rj, and also rescale D) 

^(U,Cc) = l/v^Cc + l/U-Th\ii{l/Cc) + l]h\r^{l/U) + l] + DU^, (7) 



The thermodynamical free energy F is obtained by 
minimizing with respect to £,ab and ^c- 



F{T) 



min min T{Cab,^c 

0<?a6<10<Cc<l 



(8) 



We stress again we do not claim expression (^ is a good 
detailed description of the free energy of the system, but 
only a expression (which main characteristic is given by 
Eq. (H)), that will help us to understand different sectors 
of the D-r] phase diagram of high-Tc's. 



B. Results 



We present now the results obtained by minimizing the 
free energy functional given by expression . When do- 
ing this minimization we obtain the free energy F and the 
lengths ^ab and as a function of temperature. A first 
order transition is identified as a discontinuous derivative 
of F, or equivalently a jump in the values of ^ab and ^c- 
When no jump in ^^b and is obtained, the dependence 
of £,ab and with temperature gives a clue on how the 
superconducting coherence is lost when raising tempera- 
ture. However, in this case the identification of a phase 
transition is not simple, and we can only identify temper- 
ature ranges where coherence along c- or a6-directions is 
high or low. 

Let us start with the case D — Q and rj = 1 (note 
that we are using a renormalized anisotropy parameter, 
so ?7 = 1 does not imply necessarily an isotropic system) . 
In this case Eq. (0) is symmetric in £^ab and ^c, and in 
fact the minima of are on the line ^ab = = so 
in Fig. ^ we plot the function ^(^, £,) for different val- 
ues of the temperature. We see that when T ^ the 
minimum free energy state corresponds to the maximum 
possible value of 5, i.e., the system is in the ordered state. 
When temperature increases a first order transition to a 
disordered state occurs. This is also seen from the be- 
havior of ^ as a function of temperature, as depicted in 
Fig. §(b). Thus we see that the coupling of two con- 
tinuous transition (through the entropy term in (|^)) can 
merge them onto a single one, which is first order. In 
fact, such sort of mechanism has been previously pro- 
posed to occur in other cases, such as two dimensional 
melting, in which the continuous dislocation-unbinding 
and disclination-unbinding transitions of the Kostarlitsz- 
Thouless-Halperin-Nelson- Young melting theptyOlia can 
collapse onto a first order melting transitionll3't2l 




FIG. 4. (a) Free energy functional as a function of corre- 
lation length 5 for different temperatures with D — 0, and 
7 = 2.(b) Correlation length as a function of temperature 
obtained as the minima of the curves in (a). 



If disorder is present in the system it will tend to de- 
stroy the melting transition. In Fig. |^ we show results 
as those of Fig. ^ but for a value of D = 0.7. As we 
see the jumps in and have disappeared, indicating 
that the transition is not first order. If the anisotropy 
had been chosen different from one, then the tempera- 
ture at which ^ab and take a given value would have 
been different. Although we cannot characterize from 
our simplified model a phase transition when ^ab and 
are continuous functions of temperature, it is tempting 
to say that if ^ab (Cc) drops to zero at lower temperature 
than (Cab): then we are in a sector of the phase dia- 
gram where Ti < Tp (Tp < Ti). In this way we generate 
the phase diagram depicted in Fig. ^. As indicated, it is 
qualitatively similar to the one obtained from the numer- 
ical simulation, and it gives support a posteriori to the 
proposal of an entropy of the system of the form given 
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by Eq. d). 

Some characteristics of this phase diagram can be an- 
alyzed in simple terms. For example, when 77 00, 
the system is a set of decoupled planes, and no first 
order transition is obtained for any value of D. When 
ry — !■ vortices are rigid lines and in fact effectively two- 
dimensional, so a first order transition is not obtained in 
this case either. This limiting cases give some insight on 
the form of the border between first-order and continu- 
ous transitions in Figs. |^ and ^. There is an optimum 
value of the anisotropy, at which the first order transition 
persists up to a highest value of disorder. This optimum 
value depends on the thickness of-the sample. In fact, 
as we discussed in a previous workt3 the temperature Tp 
logaritmically decreases when the thickness of the sample 
increases. This means that in Fig. ^ the border between 
the zones with Ti > Tp and Tp > Ti moves to lower val- 
ues of 77. It is thus likely that the optimum value of 77 for 
the occurrence of the first order transition also decreases 
with sample thickness . 




FIG. 5. Same as Fig. 4 but with D = 0.7. 



IV. MAGNETIC FIELD DEPENDENCE, AND 
COMPARISON WITH EXPERIMENTS 



Having discussed the D-rj phase diagram for a fixed 
magnetic field H, we turn now to the discussion of the 



dependence on H. From the numerical point of view 
the direct approach would be to do simulations at dif- 
ferent fields. However, as we discussed above, to reduce 
the magnetic field to other commensurate values would 
require exceedingly large computing time. Fortunately, 
there are arguments that suggest that a change of H can 
be mapped onto a change of D and 77. 
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FIG. 6. Disorder-anisotropy phase diagram as obtained by 
minimizing the free energy functional (see text for explana- 
tion) . 

The scaling combination between magnetic-field and 
anisotropy has been given by Chen and Teitel.tJ We gen- 
eralize here the argument to include the disorder pa- 
rameter D. Our dimensionless temperature (measured 
in units of the mean Josephson energy of the in-plane 
junctions) can be only a function of the other dimen- 
sionless parameters of the system. These are D, H, and 
77. These parameters have different dependences on the 
coherence length ^o- If we identify the discretization pa- 
rameter in the ab plane with a distance of the order of 
the coherence length ^o, then the critical current of the 
Josephson junctions along the c direction is proportional 
to ^g, so the anisotropy parameter rj behaves as ^g"^. On 
the other hand, our dimensionless magnetic field H is 
given in terms of the real external magnetic field Hq by 
the expression H = i?oCo/'^o i4>o is the flux quantum). 
For the parameter D, we note that D is proportional to 
the amplitude of the pinning potential. A vortex aver- 
ages this random function on an area ~ ^q. Considering 
the case of random (uncorrelated) pinning we find that 
D depends on ^0 as D ^ S,^^ . Since we are ignoring de- 
tails of the vortex cores, we expect the ^0 dependence to 
cancel out, and our temperature transition to be only a 
function of the ^Q-indcpendcnt quantities rj^H and D^H . 

We conclude that we can obtain the behavior of the 
system as a function of magnetic field from the results 
of Fig. ^ on lines with constant D/rj. A sketch of the 
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different possibilities is shown in Fig. |j. The general 
prediction from Fig. || is depicted in Fig. 0(a). The 
scales on the axis as well as the extent of the first order 
transition depend on the particular value of D/rj. This 
general picture has to be modified at very low fields. In 
fact, a minimum crossover fieldEJ (given essentially as the 
field at which the vortex lattice parameter matches the 
thickness of the sample) exists, below which Ti and Tp are 
essentially the same (this is because in this case there are 
so few vortices in the sample that the transition is enter- 
ally due to thermally generated vortex loops). Different 
possibilities are depicted in Fig. ^ (b), and (c). They cor- 
respond to different ranges of values of D/rj. Fig. ^(b) 
correspond to D/rj small, so Fig. ^ predicts Tp < Ti at 
high fields, and a first order zone at intermediate and 
low fields. This phase diagram corresponds to that ex- 
perimentally obtained for Bi2Sr2CaCu20,El which in fact 
has the largest value of 77, and even to the case of clean 
YBa2Cu307jH3 which has a very low D. Fig. |^(c) shows 
the expected phase diagram for samples with a high value 
oi D/rj. At low fields Ti is lower than Tp, whereas at high 
fields a cross to a zone with T, > Tp is possible. No first 
order zone is shown because a curve defined by a high 
value of D/rj in Fig. || do not pass through the first 
order zone. The low field part of this phase diagram cor- 
responds to [tie one obtained for YBa2Cu307 samples 
with defects.E3 The crossover to a case with Ti < Tp has 
not been observed, presumably because of the high fields 
needed. 




Temperature 



FIG. 7. Qualitative sketch of the H-T phase diagram for 
diferent calues of D/rj. (a) General form, (b) D/rj small, (c) 
D/rj high. 



As we indicated above, when Ti < Tp superconducting 
coherence is lost along the a6-plane at lower tempera- 
tures than along the c-axis. For Ti < T < Tp finite 
resistance within the planes and zero resistance in the c- 
dircction is observed. When Ti > Tp and in the case that 
Ti > T > Tp we potentially expect finite resistance in the 
c-direction and zero resistance within the planes. This, 
has turned to be difficult to find, both experimentallytil 
and in our simulations. The resistance seems rather to. 
go to zero at the same temperature Ti in all directions.E£l 
This is due to the fact that vortex paths crossing the 
sample for T > Tp are pinned to the a6-planes as long 
as T < Ti, preventing their movement (and thus dissi- 
pation). However, it is worth noting that some other 
experimental rrieasurements of coherence (AC magneti- 
zation) indicateE3 that in fact, c-axis coherence is lost at 
lower temperatures than in-plane coherence for the case 
of Bi2Sr2CaCu20. 

The transitions observed in our simulations are of dif- 
ferent character, and we want to discuss the point a bit 
further. The first order transition is the easiest to char- 
acterize numerically. Although we do not show all the 
results, we observed that when the resistivity has a jump 
other indicators point to a first order phase transition, 
among them the existence of hysteresis in the resistiv- 
ity curves upon heating and cooling, and the fact that 
right at the transition temperature, the energy histogram 
of the system has two peaks, indicating two coexisting 
phases with an energy barrier separating them.ESolj The 
continuous transitions are more difficult to characterize. 
The transition at Ti is not a phase transition in our 
model. In fact, it is a crossover due,,i£) thermal dep- 
pining of rather independent vorticesEiij However, in a 
real sample it may correspond to the vortex glass,transi- 
tion, depending on the strength of the disorder .lZI When 
Tp > Ti, we have previously characterized the transi- 
tion at Tp as a percolation phase transition of tlie vortex 
structure perpendicularly to the applied field.E3 In the 
thermodynamic limit for Lab {Lab ^ 00) the system does 
not have any vortex line running perpendicularly to the 
applied field for T <Tp, whereas for T > Tp these paths 
extend all over the ab plane with probability one. In Ref- 
erences [14,12] we showed numerical evidence suggesting 
that this transition is a second order phase transition and 
gave its critical exponents as found from simulations. 



V. SUMMARY AND CONCLUSIONS 

In this paper we presented numerical evidence that 
supports an anisotropy-disorder phase diagram of the 
vortex structure of high- T^s with the following char- 
acteristics: For clean samples the vortex lattice melts 
through a first order phase transition for a wide range of 
anisotropics. When disorder is included the behavior of 
the system is strongly dependent of the anisotropy. For 
low anisotropics the in-plane coherence is lost at a tem- 
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perature Ti lower than the temperature Tp at which inter- 
plane coherence is lost, and a zone of disentangled vortex 
lines is observed for Ti < T < Tp. For highly anisotropic 
samples the superconducting coherence as deduced from 
simulations of the resistivity is lost at the same temper- 
ature Ti within the planes and perpendicularly to the 
planes. However, in this case the vortex structure perco- 
lates at a temperature Tp well below Ti. In this case the 
system for Tp < T < Ti is in an 'entangled solid' phase. 
These features are also obtained from an estimation of 
the free energy of the system which is mainly based on 
a proposal for the entropy of the system (Eq. (|^)). We 
showed that the magnetic field-temperature behavior of 
the system can be deduced from results obtained from a 
fixed magnetic field provided the anisotropy and disorder 
present in the system are properly rescaled. 

Our results present in a unified way, different charac- 
teristics of the vortex structure that had been previously 
found in partial studies. The analysis is in agreement 
with a variety of experiments performed on different ma- 
terials with a broad range of parameters such as disorder, 
anisotropy and magnetic field. It could prove to be use- 
ful to find a more solid base of our proposal for the free 
energy of the system -that we showed is qualitative good- 
in order to obtain more detailed analytical results. 
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